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Abstract We give an introduction into the observational technique of integral field 
or 3D spectroscopy. We discuss advantages and drawbacks of this type of obser- 
vations and highlight a few science projects enabled by this method. In the second 
part we describe our 3D spectroscopic survey of Blue Compact Dwarf Galaxies. 
We show preliminary results from data taken with the VIMOS integral field unit 
and give an outlook on how automated spectral analysis and forthcoming instru- 
ments can provide a new view on star formation and associated processes in dwarf 
galaxies. 



1 Integral Field Spectroscopy 

Integral Field Spectroscopy (IFS) is a three-dimensional observing technique, pro- 
ducing a spectrum for every sampled element on the sky. The end result after data 
reduction is therefore typically a "datacube", with two spatial axes and the spectral 
direction. This method is therefore often called 3D Spectroscopy. 
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1.1 IFS basics 



Unlike scanning techniques, an integral field unit (IFU) splits up the telescope focal 
plane into a series of elements (typically lenslets, fibers, or slicers, or a combination 
of them) and forms a pseudoslit from which the light is dispersed to be recorded on 
a CCD. Such an exposure includes all wavelengths and full spatial sampling and ex- 
tent of the instrument, recorded at the same time and under the same (atmospheric) 
conditions. Hence, IFS is a very efficient observing scheme, especially if the objects 
of interest fit into the field of view of one pointing. This is also one of the draw- 
backs of IFS; the current instruments only cover a relatively small area on the sky, 
up to --110" (as for VIRUS-P; ||T9|). The other disadvantage is that the informa- 
tion content of each exposure is extremely high, and depending on the instrument, it 
can take a long time to properly (re-)construct the datacube from the raw data. One 
can therefore say that in general, IFS shifts the effort from the time spent to carry 
out observing programs to the time required to do proper data reduction (and data 
analysis). But recently, more and more software packages are becoming available 
for IFS data reduction. Many instrument-specific pipelines are being maintained (at 
e.g. the ESO and Gemini observatories) but there are also more general reduction 
packages for fiber-based IFU data. The newest package which has been adapted to 
reduce data from almost ten of the most popular integral field instruments is p3d[^ 

Of course, IFS would not be rising in popularity, if there were not advantages 
that by far outweigh the difficulties, at least for selected scientific applications. Ob- 
viously, it enables bidimensional spectral coverage (as opposed to traditional slit 
spectroscopy) with all the advantages that implies: no slit effects or slit losses, no 
necessity to do (pre-)imaging of the same field, a spectral resolution that is inde- 
pendent of spatial coverage, possibility to correct for atmospheric dispersion and 
to bin the data spatially to gain S/N. IFS also has potentially long wavelength cov- 
erage (as compared to the Fabry-Perot technique) and, as said above, guarantees 
homogeneous conditions for the whole dataset (as opposed to scanning techniques). 



1.2 IFS advantage: no slit effects 

The problem of slit effects [see e.g. "21 for a detailed discussion) is particularly trou- 
blesome for slit spectroscopy, which many astronomers are not aware of. If one 
observes one or more compact objects with a sUt spectrograph under seeing con- 
ditions which cause a FWHM smaller than the slit width, velocities measured for 
each object depend on their centering within the slit. This can be corrected for stellar 
fields with precisely known star positions (as done by iflTl ). But for emission line 
kinematics of compact objects such a correction is not possible, unless a narrow- 
band image with sufficient spatial resolution exists, from which the exact positions 

' Available from the project webpage http : //p3d. sourcef orge . net 
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of all relevant emission line peaks can be determined. Emission line kinematics of 
e.g. star-forming dwarf galaxies derived using long-slits is therefore to be regarded 
with caution. 

To give an example: if one observes Ha with a spectral resolution of R — 500 
at slit width of 1", the slit is effectively ^ 600 kms^' wide. At a resolution of R = 
10000 it still is - 30km s-> wide. So, if one happens to observe a compact HII 
region in 0"5 seeing with a 1" slit and R — 10000, one on average makes an error 
of ±15kms"' (for shifts of half the slit width the FWHM of the object is still 
within the slit) at each position along the slit. Rotation curves derived that way can 
look very much like disturbances caused by e.g. galaxy merging and lead to wrong 
interpretation of the galaxy at hand. 

Most kinds of integral field spectrographs are insensitive to slit effects: instru- 
ments which use lenslets are insusceptible to this problem; fiber-based systems are 
less affected as fibers scramble the spatial signal sufficiently well. Only IFUs built 
using slicers with low spatial sampling can be affected. But then one can derive the 
spatial position of most objects within the field of view, allowing to correct for the 
instrumental shifts. Of course, instruments based on the Fabry-Perot principle are 
also immune. To derive emisson line kinematics of any object one should therefore 
always choose to observe with an integral field spectrograph or a similar instrument 
that is not affected by slit effects. 



1.3 Analyses enabled by IFS 

Some kinds of analyses are only possible if one observes galaxies with IFS. Basi- 
cally, one can create two-dimensional maps of many properties, which is not possi- 
ble using slit spectroscopy. One can even create spatially resolved maps involving 
line ratios or equivalent widths (EWs) which is possible with Fabry-Perot instru- 
ments but expensive in terms of observing time. Typical IFS-based maps are those 
for extinction (Balmer decrement) and chemical abundances. Two-dimensional 
spectroscopy also allows to create color maps from which the effects of dust and 
contributions by emission line fluxes were removed; these can then be interpreted 
using photometric stellar population models without the need to build complex mod- 
els involving extinction and line emission. Kinematical analyses are of course pos- 
sible, but can be done using other observing methods. Another unique analysis en- 
abled by IFU-observations is the search for spectral signatures across a galaxy. This 
is typically done for Wolf-Rayet features in star-forming objects. 

Projects like SAURON (kinematics and stellar populations of nearby elliptical 
galaxies; |I3][T3]) and SINS (properties of z 2 galaxies; |fT5l[T4l ) would not have 
been possible without integral field spectroscopy (and a well-managed team work- 
ing on data reduction and analysis). 
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2 Studies of Blue Compact Dwarf galaxies with IFS 



Among the many recent publications involving integral field spectroscopy of dwarf 
galaxies we want to highlight those of an ongoing, comprehensive study of Blue 
Compact Dwarf (BCD) galaxies. BCDs are dwarf starburst galaxies with low metal- 
licity, often regarded as local counterparts of young galaxies observed at high- 
redshift. Although it became clear later on that the vast majority of these dwarf 
galaxies have an underlying old stellar population 19] ID, these objects are still the 
best nearby laboratories for exploring violent star formation and dwarf galaxy build 
up, including the feedback processes associated with it. Some among the many open 
questions which we try to address using IFS: What are the star- formation histories 
of BCDs? What is their dust content and distribution? Are there evolutionary con- 
nections between star-forming dwarfs and other types of dwarfs? What mechanisms 
trigger the starburst? Are interactions and mergers playing a role? 



2.1 Project description 

In our project (PI L. M. Cairos) we aim to thoroughly analyze a large sample of >40 
BCDs using integral field spectroscopy. Our targets are located in both hemispheres 
and cover a wide range in luminosity (—13.9 < Mg < —21.1 mag). We observed 
BCDs over a wide metallicity (0.05 < Zp, < 0.85) and morphology range lH. Deep 
(optical) imaging in at least two bands is available for all our targets. By now we 
have acquired data from four instruments/telescopes: INTEGRAL (coupled to the 
WYFFOS spectrograph at WHT), PMAS (3.5m Calar Alto telescope), VIRUS-? (at 
the 2.7m McDonald telescope), and VIMOS (at the ESO VLT). 

Our methodology involves a thorough analysis of the data of each object. We 
first build maps of emission line fluxes or equivalent widths for different Une types 
(hydrogen Balmer lines, forbidden oxygen or sulfur lines, etc.). These maps allow 
us to localize the HII regions and investigate the luminosity, morphology, excitation 
mechanisms, extinction pattern, and kinematics of the ionized gas component. From 
maps of the continuum we can study the morphology of the underlying stellar com- 
ponent, undisturbed by line emission, and better assess its age structure. The Balmer 
decrement enables us to carry out a spatially resolved study of the intrinsic extinc- 
tion. Further insights can be obtained from Wolf-Rayet features (based on the study 
of the emission lines in the spectral range 4650. . . 4690 A) which gives clues to the 
starburst age and the shape of the initial mass function. Metal emission lines can tell 
us about physical parameters of the gas and spatial variations of the chemical abun- 
dances. Using gas and stellar velocity fields we investigate the origin of starburst 
activity and estimate the total mass of the galaxy. Finally, we use spatially resolved 
spectral synthesis models to characterize the properties of stellar populations in the 



sample galaxies (see 3.1 



Some of our data taken with INTEGRAL and PMAS has already been published. 
An initial investigation of five galaxies with INTEGRAL was presented in lfT6ll . 
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and a paper first showing the potential of the PMAS data for the case of Mrk 1418 
appeared as Q- The PMAS data on Mrk 409 was thoroughly exploited by |5| and 
a general presentation of the objects observed so far with PMAS was published as 
H. A further description of the sub-sample of BCDs observed with the PMAS and 
VIRUS-P instruments can also be found in |[8] this volume]. All these papers show 
results that would have been impossible to acquire with other observing techniques. 



2.2 VIMOSdata 

Nine southern galaxies of our sample were so far observed with the VIMOS spectro- 
graph on the 8.2m ESO VLT. This dataset contains our deepest exposures with the 
highest spectral resolution (R ^ 2600). We used the HRblue and HRorange grisms 
to cover most of the optical wavelength range from 4300 to 7400 A. We used the 
setup with low spatial resolution (0.67" px^^) that was well matched to the seeing 
and allowed to cover ^ 27" x 27" per target. Our observing strategy involved several 
dithered exposures per galaxy and grism, to fill in bad spatial elements (spaxels) and 
to aid sky subtraction. This increased the complexity of the data reduction, since the 
combination of multiple dithered exposures into a full datacube is not supported by 
any existing processing tool. 




4500 5000 5500 6000 6500 7000 

Wavelength [Angstrom] 



Fig. 1 Spectrum of Tol 1924-416: this spectrum was integrated over the 80 central spectral ele- 
ments of our VIMOS data. It shows many faint He and metal lines with detections at good S/N. The 
lower spectrum is the same as the upper one but scaled down lOOx to show the bright emission 
lines. Grey vertical semitransparent areas indicate sky emission lines that could not be subtracted 
perfectly. 

Here, we show preliminary results of one object of our VIMOS sample, Tol 1924- 
416. Figure [T] shows an integrated spectrum of the central part of Tol 1924-416. This 
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Fig. 2 Kinematical maps of Tol 1924-416: left is the velocity field as derived from Ha, right is 
the velocity dispersion measured from the same line. 



high S/N spectrum shows many faint HeUum and metal lines that are not commonly 
seen in BCDs. Figure |2] presents maps of velocity and velocity dispersion as de- 
rived from the Ha emission line. While the emission line morphology is similar 
to the broad-band appearance, the velocity structure is very complex. While local 
starburst-driven winds might contribute to this structure, it is tempting to take it as 
a sign of merging, similar to what previous investigations of this galaxy concluded 
from optical broad- and narrow-band observations ||T8][T2|. Derivation of stellar ve- 
locities in this galaxy, to confirm or disprove this view, are difficult as only very few 
spaxels show absorption. 



3 Conclusions & Outlook 

In this article we described the technique of integral field spectroscopy that provides 
an unbiased and unrestricted view on dwarf galaxies: from the same data, taken 
under nearly identical conditions, one can derive kinematics (gaseous and possibly 
stellar) and physical properties of the warm interstellar medium (e.g. metallicity, 
density), and study the stellar populations in two dimensions. We briefly outlined 
a project investigating the structure and evolution of BCDs with IFS and showed 
preliminary results for one galaxy from our sample. 

But most of the findings enabled by 3D spectroscopy are probably still to come. 
Next, we briefly discuss two topics that will help to further exploit existing data and 
to produce major new results, in the research area of dwarf galaxies and beyond. 
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Fig. 3 Automated stellar continuum analysis of VIRUS-P data for the BCD IlIZw 102: the left 
and right panels show, respectively, the H/3 map of the galaxy prior to and after correction for 
underlying stellar absorption. 



3.1 Automated spectral analysis 

Most IFS-based investigations so far analyze the integrated spectrum of multiple 
regions of interest, or perform Gaussian emission line fitting to individual spaxels, 
leaving unexploited the information provided by stellar absorption lines. As deeper 
data gets available, the latter will become more accessible and can be used, espe- 
cially for kinematical studies. 

Another approach is to use spectral synthesis and analysis tools to model the 
stellar continuum spaxel by spaxel. This yields a twofold benefit, as it allows for 
insights into the stellar populations and for the correction for underlying stellar ab- 
sorption, thereby an improvement on the accuracy of emission line measurements. 
Especially Balmer line fluxes usually have to be corrected for underlying stellar ab- 
sorption. The subtraction of the best-fitting stellar model offers an effective means 
for accomplishing this task. 

As an example, we show in Fig. |3] preliminary results from the analysis of 
our VIRUS-P data for the BCD IIIZw 102 (Cairos et al., in prep.). We used the 
STARLIGHT code ifToHTTI to generate a stellar continuum fit at each position. Fig- 
ure[3]shows the H)3 flux before and after correction for underlying stellar absorption 
using the best-fitting stellar model for each spaxel. It can be seen that the artifical 
depression in the Hj3 flux distribution in southeastern and northwestern direction 
could be efficiently rectified this way (Papaderos et al., in prep.). 



8 Peter M. Weilbacher, Luz Marina Cairos, Nicola Caon, Polychronis Papaderos 

3.2 New instruments 



Several major new IFSs are scheduled to come on-line in the next years. At the 
ESO VLT, all new instruments will be based on the integral field principle or will 
at least contain an ancillary IFU-mode. The major IFSs to become available at ESO 
are KMOS, a near-infrared spectrograph with 24 small IFUs, and MUSE, an optical 
instrument. The VIRUS instrument in development for the HET at McDonald ob- 
servatory will provide simultanous (but not contiguous) spectral coverage of an area 
of 30 □' using a setup of up to ^ 150 fiber bundles. 

Of these, the Multi Unit Spectroscopic Explorer (MUSE; [IJ) will probably be 
the most powerful instrument to further drive ahead the field of nearby star-forming 
dwarf galaxies. It is currently being built as one of the 2nd generation instruments 
for the ESO VLT and will see first fight in 2012. Its I'x 1' field of view at 0"2 
sampling, long wavelength range of (4650) 4800. . .9300 A, relatively high spectral 
resolution {R ^ 3000), and high throughput should make it an ideal multi-purpose 
instrument. After first light, it will be enhanced with an adaptive optics module to in 
the end allow near-diffraction limited spectroscopy in a 7"5x7"5 field with a 0"025 
sampling. 
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